SIMU=0;CD_REAL_5_0=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_0[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_0=c(CD_REAL_5_0,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_0,'CD_REAL_5_0.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_1=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_1[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_1=c(CD_REAL_5_1,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_1,'CD_REAL_5_1.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_2=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_2[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_2=c(CD_REAL_5_2,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_2,'CD_REAL_5_2.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_3=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_3[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_3=c(CD_REAL_5_3,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_3,'CD_REAL_5_3.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_4=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_4[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_4=c(CD_REAL_5_4,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_4,'CD_REAL_5_4.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_5=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_5[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_5=c(CD_REAL_5_5,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_5,'CD_REAL_5_5.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_6=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_6[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_6=c(CD_REAL_5_6,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_6,'CD_REAL_5_6.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_7=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_7[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_7=c(CD_REAL_5_7,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_7,'CD_REAL_5_7.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_8=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_8[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_8=c(CD_REAL_5_8,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_8,'CD_REAL_5_8.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_9=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_9[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_9=c(CD_REAL_5_9,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_9,'CD_REAL_5_9.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_10=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_10[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_10=c(CD_REAL_5_10,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_10,'CD_REAL_5_10.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_15=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_15[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_15=c(CD_REAL_5_15,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_15,'CD_REAL_5_15.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_20=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_20[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_20=c(CD_REAL_5_20,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_20,'CD_REAL_5_20.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_25=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_25[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_25=c(CD_REAL_5_25,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_25,'CD_REAL_5_25.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_30=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_30[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_30=c(CD_REAL_5_30,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_30,'CD_REAL_5_30.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_35=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_35[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_35=c(CD_REAL_5_35,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_35,'CD_REAL_5_35.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_40=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_40[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_40=c(CD_REAL_5_40,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_40,'CD_REAL_5_40.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_45=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_45[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_45=c(CD_REAL_5_45,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_45,'CD_REAL_5_45.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_50=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_50[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_50=c(CD_REAL_5_50,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_50,'CD_REAL_5_50.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_55=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_55[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_55=c(CD_REAL_5_55,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_55,'CD_REAL_5_55.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_60=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_60[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_60=c(CD_REAL_5_60,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_60,'CD_REAL_5_60.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_65=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_65[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_65=c(CD_REAL_5_65,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_65,'CD_REAL_5_65.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_70=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_70[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_70=c(CD_REAL_5_70,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_70,'CD_REAL_5_70.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_75=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_75[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_75=c(CD_REAL_5_75,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_75,'CD_REAL_5_75.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_80=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_80[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_80=c(CD_REAL_5_80,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_80,'CD_REAL_5_80.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_85=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_85[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_85=c(CD_REAL_5_85,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_85,'CD_REAL_5_85.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_90=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_90[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_90=c(CD_REAL_5_90,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_90,'CD_REAL_5_90.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_95=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_95[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_95=c(CD_REAL_5_95,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_95,'CD_REAL_5_95.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_100=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_100[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_100=c(CD_REAL_5_100,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_100,'CD_REAL_5_100.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_110=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_110[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_110=c(CD_REAL_5_110,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_110,'CD_REAL_5_110.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_120=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_120[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_120=c(CD_REAL_5_120,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_120,'CD_REAL_5_120.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_130=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_130[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_130=c(CD_REAL_5_130,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_130,'CD_REAL_5_130.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_140=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_140[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_140=c(CD_REAL_5_140,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_140,'CD_REAL_5_140.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_150=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_150[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_150=c(CD_REAL_5_150,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_150,'CD_REAL_5_150.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_160=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_160[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_160=c(CD_REAL_5_160,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_160,'CD_REAL_5_160.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_170=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_170[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_170=c(CD_REAL_5_170,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_170,'CD_REAL_5_170.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_180=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_180[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_180=c(CD_REAL_5_180,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_180,'CD_REAL_5_180.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_190=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_190[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_190=c(CD_REAL_5_190,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_190,'CD_REAL_5_190.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_200=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_200[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_200=c(CD_REAL_5_200,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_200,'CD_REAL_5_200.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_210=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_210[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_210=c(CD_REAL_5_210,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_210,'CD_REAL_5_210.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_220=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_220[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_220=c(CD_REAL_5_220,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_220,'CD_REAL_5_220.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_230=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_230[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_230=c(CD_REAL_5_230,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_230,'CD_REAL_5_230.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_240=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_240[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_240=c(CD_REAL_5_240,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_240,'CD_REAL_5_240.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_250=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_250[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_250=c(CD_REAL_5_250,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_250,'CD_REAL_5_250.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_260=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_260[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_260=c(CD_REAL_5_260,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_260,'CD_REAL_5_260.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_270=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_270[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_270=c(CD_REAL_5_270,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_270,'CD_REAL_5_270.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_280=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_280[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_280=c(CD_REAL_5_280,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_280,'CD_REAL_5_280.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_290=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_290[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_290=c(CD_REAL_5_290,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_290,'CD_REAL_5_290.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_300=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_300[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_300=c(CD_REAL_5_300,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_300,'CD_REAL_5_300.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_310=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_310[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_310=c(CD_REAL_5_310,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_310,'CD_REAL_5_310.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_320=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_320[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_320=c(CD_REAL_5_320,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_320,'CD_REAL_5_320.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_330=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_330[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_330=c(CD_REAL_5_330,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_330,'CD_REAL_5_330.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_340=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_340[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_340=c(CD_REAL_5_340,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_340,'CD_REAL_5_340.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_350=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_350[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_350=c(CD_REAL_5_350,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_350,'CD_REAL_5_350.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_360=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_360[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_360=c(CD_REAL_5_360,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_360,'CD_REAL_5_360.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_370=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_370[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_370=c(CD_REAL_5_370,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_370,'CD_REAL_5_370.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_380=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_380[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_380=c(CD_REAL_5_380,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_380,'CD_REAL_5_380.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_390=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_390[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_390=c(CD_REAL_5_390,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_390,'CD_REAL_5_390.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_400=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_400[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_400=c(CD_REAL_5_400,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_400,'CD_REAL_5_400.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_410=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_410[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_410=c(CD_REAL_5_410,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_410,'CD_REAL_5_410.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_420=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_420[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_420=c(CD_REAL_5_420,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_420,'CD_REAL_5_420.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_430=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_430[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_430=c(CD_REAL_5_430,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_430,'CD_REAL_5_430.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_440=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_440[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_440=c(CD_REAL_5_440,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_440,'CD_REAL_5_440.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_450=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_450[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_450=c(CD_REAL_5_450,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_450,'CD_REAL_5_450.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_460=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_460[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_460=c(CD_REAL_5_460,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_460,'CD_REAL_5_460.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_470=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_470[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_470=c(CD_REAL_5_470,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_470,'CD_REAL_5_470.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_480=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_480[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_480=c(CD_REAL_5_480,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_480,'CD_REAL_5_480.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_490=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_490[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_490=c(CD_REAL_5_490,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_490,'CD_REAL_5_490.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5_500=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_500[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5_500=c(CD_REAL_5_500,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5_500,'CD_REAL_5_500.txt',row.names=FALSE,col.names=FALSE)

